System size dependence of elliptic flows in relativistic heavy-ion collisions 
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The elliptic flows in both Cu+Cu and Au+Au collisions at the Relativistic Heavy Ion Collider are 
studied in a multi-phase transport model. For both collisions at same reduced impact parameter and 
minimum bias collisions, the elliptic flow of partons in Cu+Cu collisions is about a factor of three 
smaller than that in Au+Au collisions at same energy. The reduction factor is similar to the ratio 
of the sizes of the two colliding systems and is also related to the combined effects of initial energy 
density and spatial elliptic deformation in the two reactions. A similar system size dependence is 
also seen in the elliptic flow of charged hadrons from minimum bias collisions. 



PACS numbers: 25.75.Ld, 24.10.Lx 

I. INTRODUCTION 

There have been extensive studies on the azimuthal 
anisotropy of hadron momentum distributions in the 
transverse plane perpendicular to the beam direction, 
particularly the elliptic flow V2 , in heavy-ion collisions at 
various energies EI- The hadron transverse momentum 
anisotropy is generated by the pressure anisotropy in the 
initial compressed matter formed in non-central heavy 
ion collisions HH and thus depends on the initial geome- 
try and energy density of a collision as well as the proper- 
ties of produced matter during the early stage of the col- 
lision. For heavy-ion collisions at the Relativistic Heavy 
Ion Collider (RHIC).it has been shown that not only the 
larger elliptic flow 0, H H 0, II E| but also the smaller 
higher-order anisotropic flows Halll El El El El 13 
are sensitive to the properties of initial dense matter. To 
study the dependence of the hadron anisotropic flow on 
the initial geometry and energy density in a heavy ion 
collision, one can vary the impact parameter of the colli- 
sion or the atomic number of colliding nuclei. Although 
there were already extensive studies on the dependence 
of elliptic flow on the impact parameter of heavy ion col- 
lisions, its dependence on the size of colliding nuclei has 
only been started recently Qjj . 

In the present work, we use a multi-phase transport 
(AMPT) model, that includes both initial partonic and 
final hadronic interactions [3 El, to study the elliptic 
flow in both Cu+Cu and Au+Au collisions at s/s — 200 
AGeV at RHIC. In particular, we use the version with 
string melting, i.e., allowing hadrons that are expected to 
be formed from initial strings to convert to their valence 
quarks and antiquarks |2Ct l2ll |22j , which has been shown 
to be able to explain the measured pr dependence of «2 
and Vi of mid-rapidity charged hadrons with a parton 
scattering cross section of about 10 mb. We find that 
the elliptic flow scales approximately with the size of the 
colliding system and is also proportional to the product of 
the spatial elliptic deformation and energy density during 
the initial stage of a collision. 

This paper is organized as follows. In Sec. [HI the 



AMPT model is briefly reviewed. Results for the parton 
elliptic flow and spatial elliptic deformation as well as 
their system size dependence are shown in Sec. 11111 while 
the system size dependence of the elliptic flow of charged 
hadrons is shown in Sec. IIVI Finally, a brief summary is 
given in Sec. IVl 

II. THE AMPT MODEL 

The AMPT model E1E100H1 is a hybrid model 
that uses minijet partons from hard processes and strings 
from soft processes in the Heavy Ion Jet Interaction Gen- 
erator (HIJING) model 2@J as the initial conditions for 
modeling heavy ion collisions at ultra-relativistic ener- 
gies. Since the initial energy density in Au+Au collisions 
at RHIC is much larger than the critical energy density at 
which the hadronic matter to quark-gluon plasma transi- 
tion would occur |23ll27|. we use the version which allows 
the melting of initial excited strings into partons [2(J • In 
this string melting scenario, hadrons (mostly pions), that 
would have been produced from string fragmentation, are 
converted instead to valence quarks and/or antiquarks 
with current quark masses. Interactions among these 
partons are described by Zhang's parton cascade (ZPC) 
model |28|. At present, this model includes only parton- 
parton elastic scatterings with an in-medium cross sec- 
tion given by: 

dt ~ 2 \ s J (i-M 2 ) 2 ' 1 J 

where the strong coupling constant a s is taken to be 0.47, 
and s and t are usual Mandelstam variables. The effec- 
tive screening mass (X depends on the temperature and 
density of the partonic matter but is taken as a parameter 
in ZPC for fixing the magnitude and angular distribution 
of parton scattering cross section. Since there are no in- 
elastic scatterings, only quarks and antiquarks from the 
melted strings are present in the partonic matter. The 
transition from the partonic matter to the hadronic mat- 
ter is achieved using a simple coalescence model, which 
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combines two nearest quark and antiquark into mesons 
and three nearest quarks or antiquarks into baryons or 
anti-baryons that are close to the invariant mass of these 
partons. The present coalescence model is thus some- 
what different from the ones recently used extensively 
Hil EH f° r studying hadron production at inter- 
mediate transverse momenta. The final-state hadronic 
scatterings are then modeled by a relativistic transport 
(ART) model j3^|. Using parton scattering cross sections 
of 6-10 mb, the AMPT model with string melting is able 
to reproduce both the centrality and transverse momen- 
tum (below 2 GeV/c) dependence of the elliptic flow [2(j 
and pion interferometry plj measured in Au+Au colli- 
sions at Vs = 130 AGeV at RHIC El- It has also 
been used for studying the kaon interferometry |36j | and 
the elliptic flow of charmed mesons [31} in these collisions. 
We note that the above cross sections are significantly 
smaller than that needed to reproduce the parton ellip- 
tic flow from the hydrodynamic model [3^. The result- 
ing hadron elliptic flows in the AMPT model with string 
melting are, however, amplified by modeling hadroniza- 
tion via quark coalescence ^S^, leading to a satisfactory 
reproduction of experimental data. 



III. ELLIPTIC FLOW AND SPATIAL 
ANISOTROPY OF PARTONS 

The anisotropic flows v n of particles are the Fourier 
coefficients in the decomposition of their transverse mo- 
mentum spectra in the azimuthal angle (f> with respect to 
the reaction plane 39], i.e., 
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[l + J22v n (p T ,y)cos(nct>)} (2) 



n=l 



Because of the symmetry —0 in collision geome- 

try, no sine terms appear in above expansion. For parti- 
cles at midrapidity in collisions with equal mass nuclei, 
anisotropic flows of odd orders vanish as a result of the 
additional symmetry «-> (f> + tt. Anisotropic flows gen- 
erally depend on particle transverse momentum and ra- 
pidity, and for a given rapidity the anisotropic flows at 
transverse momentum px can be evaluated according to 



v n {pr) = (cos(n0)) 



(3) 



where (• • •) denotes average over the azimuthal distri- 
bution of particles with transverse momentum px- The 
elliptic flow V2 can further be expressed in terms of single- 
particle averages: 



V2{pt) 




(4) 



where p x and p v are, respectively, projections of the par- 
ticle momentum in and perpendicular to the reaction 
plane. 



Since the AMPT model also provides information on 
the spatial anisotropy of colliding matter, which is re- 
sponsible for generating the momentum anisotropic flows, 
it is of interest to introduce the spatial anisotropic coeffi- 
cient s n by expressions similar to those for the anisotropic 
flows v n but in terms of the spatial distributions of par- 
ticles in the transverse plane. In particular, the spatial 
elliptic deformation S2 or eccentricity can be obtained 
from 
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(5) 



where x and y are, respectively, projections of the particle 
coordinate in and perpendicular to the reaction plane. 
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FIG. 1: (Color online) Time evolution of spatial elliptic de- 
formation S2 and momentum elliptic flow v% of partons at 
midrapidity in Cu+Cu collisions at b = 5.5 fm (open squares) 
and Au + Au collisions at b = 8 fm (solid squares) with the 
same energy per nucleon y^s = 200 AGeV for parton scatter- 
ing cross sections tj p = 3 mb (left panels) and a v = 10 mb 
(right panels). 

Using the AMPT model in the string melting scenario 
with parton scattering cross sections a v — 3 and 10 mb, 
we first study the time evolution of S2 and V2 of partons 
at midrapidity in Cu+Cu collisions at b = 5.5 fm and 
Au+Au collisions at b = 8 fm with the same energy per 
nucleon, i.e., y/s = 200 AGeV, and the results are shown 
in Fig. ^ The two impact parameters, which are similar 
to the average values in minimum bias collisions, cor- 
respond to the same reduced impact parameter &/& max , 
where & max is the sum of the radii of colliding nuclei. 
They would lead to same S2 if the two nuclei were sharp 
spheres. Because of diffused surfaces, a much large S2 is, 
however, produced in Au+Au collisions than in Cu+Cu 
collisions as shown in Figs. H^a) and^b). For both colli- 
sions, the spatial elliptic deformation s 2 is initially large 
and decreases with time. It reaches the saturation value 
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at times which depend on both the size of the reaction 
system and the parton scattering cross section. Specifi- 
cally, the S2 in Cu+Cu collisions has an earlier saturation 
time, implying that the fireball (or quark gluon plasma) 
formed in a smaller reaction system has a shorter life- 
time. This behavior is understandable since the larger 
reaction system has a larger initial parton number den- 
sity in the transverse plane and thus takes longer time 
to freeze out. The nonzero spatial elliptic deformation 
S2 indicates that the parton spatial distribution is non- 
spherical at freeze out, with the larger parton scattering 
cross section leading to a smaller spatial anisotropy. 

Time evolution of the elliptic flow vi of partons is 
shown in Figs. Wi c ) an d E^d) for the two parton scat- 
tering cross sections of 3 and 10 mb, respectively. The 
elliptic flow V2 is seen to saturate earlier in the colli- 
sions, and the saturation time depends on the size of the 
reaction system, i.e., at about 5 fm/c for Cu+Cu colli- 
sions while about 7 fm/c for Au+Au collisions. This is 
similar for both parton scattering cross sections, except 
that the larger one leads to a larger elliptic flow. The 
shorter saturation time of in lighter reaction system 
is consistent with the earlier saturation time of S2 seen 
in Figs, da) an d Wfo) for that system. An interesting 
result predicted by the AMPT model is that the parton 
V2 in Cu+Cu collisions is significantly smaller than that 
in Au+Au collisions. For instance, the final values of V2 
of partons are about 1.9% and 2.9% , respectively, for 
<7 p = 3 mb and 10 mb in Cu+Cu collisions, and they in- 
crease to about 5.4% and 7.5% , respectively, for <j p — 3 
mb and 10 mb in Au+Au collisions. For both cross sec- 
tions, the strength of parton V2 is enhanced by a factor 
of about 3 from Cu+Cu to Au+Au collisions at same 
reduced impact parameter, and is roughly similar to the 
ratio of the sizes of Au and Cu nuclei. 

The magnitude of elliptic flow is affected not only by 
the spatial elliptic deformation but also by the energy 
density during the initial stage of collisions. In Fig. 
we show the time evolution of the energy densities of 
midrapidity partons (upper panel) and their ratio (lower 
panel) in Cu+Cu collisions at b = 5.5 fm (dashed line) 
and Au+Au collisions at b = 8 fm (solid line ) fo r parton 
scattering cross section of 10 mb. As in Ref. [23|, the en- 
ergy density is calculated for partons in the central cell 
that has a transverse radius of 1 fm and a longitudinal 
dimension of 5% of the time t with the time coordinate 
starting when the two nuclei fully overlap in the longitu- 
dinal direction. It is seen that the initial energy density 
in Cu+Cu collisions is about 2/3 of that in Au+Au col- 
lisions. Combining this effect with that due to initial 
spatial elliptic deformation, which has a ratio of about 
1/2 for the two colliding systems, gives a reduction fac- 
tor of about 1/3, which is comparable to the ratio of the 
sizes of the two colliding systems. The dependence on 
the system size seen in the parton elliptic flows from the 
AMPT model is thus related to the combined effects of 
initial energy density and spatial elliptic deformation. 

The same system size dependence is also seen in the 
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FIG. 2: (Color online) Time evolution of energy densities 
of midrapidity partons (upper panel) and their ratio (lower 
panel) in Cu+Cu collisions at 6 = 5.5 fm (dashed line) and 
Au + Au collisions at b = 8 fm (solid line) with the same 
energy per nucleon ^Js — 200 AGeV for parton scattering 
cross sections a v = 10 mb. 
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FIG. 3: (Color online) Transverse momentum pr dependence 
of «2 of midrapidity partons from minimum bias events in 
Au+Au (solid squares) and Cu+Cu (open squares) collisions 
at = 200 AGeV for parton scattering cross sections of 
o p — 3 mb (left panel) and a v = 10 mb (right panel). The 
solid line is 0.3 (0.35) times vi of Au+Au for a v = 3 (10) mb. 



parton differential elliptic flow as a function of parton 
transverse momentum px for minimum bias collisions. 
This is shown in Figs. E^a) and[3b) for midrapidity par- 
tons from minimum bias events in Au+Au and Cu+Cu 
collisions at ^/s — 200 AGeV for parton scattering cross 
sections of 3 mb and 10 mb, respectively. As shown by 
the solid lines in Figs. C^a) and|3fb), the parton differen- 
tial V2 in Cu+Cu collisions is about 0.3 (0.35) times that 
in Au+Au collisions for parton scattering cross section 
of 3 (10) mb. 
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We note that a detailed comparison between the par- 
ton elliptic flows in the two colliding systems further in- 
dicates that the parton v-i in the lighter reaction system 
is more sensitive to the parton cross sections. For exam- 
ple, the parton v-i is enhanced by a factor of about 1.53 
in Cu+Cu collisions but is enhanced by a factor of about 
1.39 in Au+Au collisions, when the parton cross section 
is increased by a factor of about 3. 



IV. SYSTEM SIZE DEPENDENCE OF THE 
CHARGED HADRON ELLIPTIC FLOW 
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FIG. 4: (Color online) Same as Fig. [3] but for charged 
hadrons. The solid line is 0.35 (0.4) times vi of Au+Au for 
cr p = 3 (10) mb. 

The anisotropic flows of partons are transferred to 
those of hadrons when the latter are formed from the 
coalescence of quarks and/or antiquarks. Although scat- 
terings among hadrons are included in the AMPT model, 
they do not affect much the hadron anisotropy flows as 
a result of the small spatial anisotropy and pressure in 
the hadronic matter |20j . In Fig. we show the fi- 
nal i>2 of charged hadrons in the pseudorapidity range 
< 1.2 in minimum bias Au+Au and Cu+Cu collisions 
at yfs — 200 AGeV as functions of the transverse momen- 
tum pt for parton scattering cross sections o v — 3 and 10 
mb. It is seen that the charged hadron V2 in the lighter 
reaction system Cu+Cu is significantly smaller than that 
in the heavier reaction system Au+Au. As shown by the 
solid lines in Figs. Ela) and^Jb), the charged hadron 
t>2 in Cu+Cu collisions is about 0.35 (0.4) times that in 
Au+Au collisions for parton scattering cross section of 3 
(10) mb. The elliptic flow of charged hadrons thus shows 
a similar dependence on the size of the colliding system 
as that of partons. 

A similar dependence on the size of collision system 
is also seen in the transverse momentum distribution 
of charged hadrons as shown in Fig. [S] For charged 
hadrons in the pseudorapidity range I77I < 1.2 in mini- 
mum bias Au+Au and Cu+Cu collisions at yfs = 200 
AGeV for parton scattering cross sections a v = 3 and 
10 mb, the transverse momentum spectrum of charged 



10 

a- 10' 

> 

O 10 

™o 

H-10 1 
CL 
TJ 

Z 10° 

sTio"' 

10" 2 
10"' 



CM 



: (a) o =3 mb 


(b)o = 10mb ' 


* AMPT with string melting 




| ■ Mini, bias at s 1,2 =200 AGeV 




a^^x Charged hadrons at midrapidity 












: — ■— Au+Au ^^NT^* 


■ Au+Au >v X " 


] —a— Cu+Cu ^sT^l 


■ -o- Cu+Cu >v i 


: Au+Au (*0.32) > 


Au+Au (*0.32) > 



0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 2.5 

P T (GeV/c) 

FIG. 5: (Color online) Transverse momentum distribution of 
charged hadrons in the pseudorapidity range \rj\ < 1.2 from 
minimum bias events in Au+Au (solid squares) and Cu+Cu 
(open squares) collisions at yfs = 200 AGeV for parton scat- 
tering cross sections of a v — 3 mb (left panel) and a v = 10 
mb (right panel). Solid lines are the results from Au+Au 
collisions that are scaled by 0.32. 



hadrons in Cu+Cu collisions is again about 0.32 times 
that in Au+Au collisions as shown by the solid lines in 
Figs. I3a) and[5fb). This implies that the total trans- 
verse energy (uit) of hadrons at midrapidity follows the 
same dependence on the system size. Since the initial 
energy density in heavy ion collisions can be estimated 
via e ~ (itit) /(tottR 2 ) |40| . where To is the formation 
or thermalization time and R is the radius of the collid- 
ing nuclei, the ratio of its values in Cu+Cu and Au+Au 
collisions is thus 



£Cu+Cu 
£Au+Au 



(ra T ) Cu+C u 

(TOT) Au+Au 



0.68. 



(6) 



This estimated value is consistent with that from the 
AMPT model shown in Fig. 



V. SUMMARY 

Using the AMPT model with string melting, we have 
studied elliptic flows in Cu+Cu and Au+Au collisions 
at RHIC. We find that for both collisions at same re- 
duced impact parameter and minimum bias collisions, 
the elliptic flow of partons in the lighter Cu+Cu colli- 
sions is about a factor of 3 smaller than that in the heav- 
ier Au+Au collisions at same energy, similar to the ratio 
of the system sizes. By examining the energy density and 
spatial elliptic deformation during the initial stage of a 
collision, we further find that their combined effect is sim- 
ilar to the size effect seen in the elliptic flows in the two 
colliding systems. A similar ratio is predicted for the el- 
liptic flows of charged hadrons in minimum bias collisions 
of the two systems. Recent preliminary experimental re- 
sults from the PHOBOS collaboration [17| have shown 
that if the charge hadron elliptic flow is scaled by the 
initial eccentricity or spatial elliptic deformation of par- 
ticipant nucleons, its magnitude becomes qualitatively 
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similar in both Cu+Cu and Au+Au collisions. How our 
results are related to this observation is of great interest 
and will be addressed in a subsequent study. 
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